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Localization transitions as a function of temperature require a many-body mobility edge in energy, 
separating localized from ergodic states. We argue that this scenario is inconsistent because local 
fluctuations into the ergodic phase within the supposedly localized phase can serve as mobile bubbles 
that induce global delocalization. Such fluctuations inevitably appear with a low but hnite density 
anywhere in any typical state. We conclude that the only possibility for many-body localization 
to occur are lattice models that are localized at all energies. Building on a close analogy with a 
model of assisted two-particle hopping, where interactions induce delocalization, we argue why hot 
bubbles are mobile and do not localize upon diluting their energy. Numerical tests of our scenario 
show that previously reported mobility edges cannot be distinguished from hnite-size effects. 

PACS numbers: 05.30.Rt, 72.15.Rn, 72.20.Ee 


I. INTRODUCTION 

It is now almost mathematically proven that many- 
body localization, i.e., the absence of long-range trans¬ 
port in a thermodynamic many-body system, occurs in 
certain one-dimensional quantum lattice models at any 
energy density if sufficiently strong quenched disorder is 
present [l|. In this case, many-body localization (MBL) 
comes along with a complete set of conserved quasi¬ 
local quantities However, it remains less clear 

whether the originally predicted localization transition 
at finite temperature [3, 0 exists as a genuine dynam¬ 
ical phase transition defining a sharp many-body mo¬ 
bility edge in energy density. Even though several nu¬ 
merical investigations in small one-dimensional (ID) sys¬ 
tems have reported such mobility edges , studies 
in larger systems did not find similar evidence [Mini 
and, moreover, linked-cluster analysis [a of the numer¬ 
ical data hint that the extent of the localized phase has 
been vastly overestimated. Furthermore, recent theoret¬ 
ical considerations [3, IT3I4T0 have raised doubts about 
non-perturbative effects which might reduce the putative 
transition to a crossover. A related open issue concerns 
the many-body analog of Mott’s argument, which forbids 
the coexistence of localized and delocalized states at the 
same energy in single-particle problems. 

In this paper, we address these issues, which are fun¬ 
damental for a complete understanding of localization, 
equilibration, and transport in closed many-body quan¬ 
tum systems. We argue that for systems with short- 
range interactions, many-body mobility edges cannot ex¬ 
ist, thus ruling out sharp transitions from a conducting 
to a completely insulating phase as a function of temper¬ 


ature. These considerations also imply a strong many- 
body version of Mott’s argument, which rules out the 
coexistence of localized and delocalized states, even at 
extensively different energies. A simple generalization of 
our considerations rules out mobility edges as a function 
of any extensive thermodynamic parameter, such as par¬ 
ticle number or magnetization. 


Our paper is organized as follows. We introduce all 
our arguments at a non-technical level in Sec. im and 
argue that local hot thermal spots, dubbed bubbles, con¬ 
stitute a mechanism for global delocalization. Section Hill 
contains a more detailed presentation of the argument, 
while Sec. m is devoted to the analysis and discussion 
of potential caveats. Our numerical results are presented 
in Sec. El By a careful analysis of the model considered 
by Kjall et al.^, we show that currently available sys¬ 
tem sizes are too small to host a truly thermal bubble, 
and, hence, that existing numerical data do not contra¬ 
dict our theory. To conclude, in Sec. [Vll we discuss the 
physical consequences of our analysis, and point out, in 
which physical systems it allows a genuine MBL phase to 
exist. More detailed information and discussions are rel¬ 
egated to three appendixes: Appendix [Al presents numer¬ 
ical results demonstrating delocalization via rare events 
in a two-particle model. In Appendixes [B] and o we 
discuss why the idea of bubbles as rare fluctuations of 
energy or density neither apply to many-body systems 
with disorder which localizes the full spectrum, nor to 
single-particle problems. 
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II. SUMMARY OF THE ARGUMENTS 

In this section, we first discuss a simple model of as¬ 
sisted two-particle hopping, which illustrates several im¬ 
portant features that this problem has in common with 
the rare events that induce delocalization in many-body 
systems (see Appendix for numerical results). For 
this model, we show how rare local fluctuations induce 
hybridization among putatively localized states. Then, 
we turn to our main topic, general many-body systems. 
We explain in non-technical terms how those rare events 
wash out mobility edges whenever there is an ergodic 
state at some finite energy or particle density (or, more 
generally, if there are ergodic states in any region of the 
parameter space for the extensive thermodynamic quan¬ 
tities). Finally, we discuss some potential caveats, and 
argue why they are benign. 


A. Assisted hopping model 

Consider particles in a hypercubic lattice of linear size 
L, hopping with amplitude ti between nearest-neighbor 
sites, and subject to a disorder potential i.i.d. uni¬ 
formly in [—W^W]. A particle on site x interacts with 
others by inducing an assisted hopping of strength t 2 
along the diagonals of plaquettes that x belongs to 

H =-ti^{clcy + R.c.) + '^e^na; (1) 

{x,y) ® 

-t2j2 +H.C.), 

X s,s'=±l l<a<l3<d 

where are lattice unit vectors. This may describe 

the effect of a lattice distortion brought about by the 
presence of the first particle. This model is illustrated in 

Fig. m 
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a composite light particle and overcome Anderson local¬ 
ization. This effect is related to the interaction-induced 
increase of the localization length in sufficiently weakly 
localized systems [illii. A single-particle analog of the 
phenomenon is the solvable case of two coupled Bethe 
lattices El- The delocalization in o seems natural, 
since all configurations of two particles at distance one 
are strongly resonant with each other. They thus form 
a percolating, delocalized resonant subgraph in configu¬ 
ration space, which supports delocalized wave functions 
with inverse participation ratios that vanish as the in¬ 
verse volume. We have numerically confirmed this delo¬ 
calization effect [see Appendix (see also Refs. HiHii 
for similar observations)]. In a system of only two parti¬ 
cles, the eigenstates come in two kinds: the overwhelming 
number of states is strongly concentrated on a configu¬ 
ration with two distant, immobile particles. Only a van¬ 
ishing fraction of order [log{L)]^ /L of all two-particle 
eigenstates are delocalized as dynamically bound, mo¬ 
bile pairs. As we will discuss in more detail in Sec. uni 
it would be misleading to think that two particles that 
start off together will always lose sight of each other after 
some time and localize far from each other. Instead, we 
will argue that the delocalization channel via pair config¬ 
urations, and thus delocalized eigenfunctions, are robust. 

In this two-particle model, localized and delocalized 
states coexist at the same energy. In contrast to generic 
single particle problems, this is possible here because 
the matrix elements that couple the two kinds of states 
through a random perturbation of the Hamiltonian are 
typically exponentially small in the system size. They are 
thus negligible as compared to the relevant level spacings 
and hence do not hybridize the two types of states. 

Let us now discuss how a finite density of particles 
modifies the situation. In the thermodynamic limit, there 
is a finite density of close pairs in typical configurations. 
These pairs diffuse through the sample. Initially well- 
isolated and localized particles scatter inelastically off 
these pairs and thus move as well, leading to complete 
delocalization. Even in exponentially rare configurations 
where initially all particles are far from each other, par¬ 
ticles eventually tunnel together and decay into the con¬ 
tinuum of diffusive pair states. We thus do not expect 
any localized eigenstates to survive at finite density. 


FIG. 1: Hopping processes allowed by the Hamiltonian of 
Eq. ([1]) in d = 2. Left panel: single-particle hopping. Right 
panel: assisted hopping for the left particle due to the 
presence of the right particle. 

Let us first focus on the particular case where there are 
only two particles in the system. We consider parameters 
ti <C W, for which the single-particle problem is local¬ 
ized in the whole spectrum. Eor ^2 ^ W, the two-particle 
problem has several interesting features: In dimensions 
d > 2, the assisted hopping term induces a delocaliza¬ 
tion of close pairs which will move together diffusively as 


B. Many-body systems: 

Delocalization from rare bubbles 

The argument of Basko et al @ for a localization tran¬ 
sition as a function of temperature, i.e., a many-body 
mobility edge, builds on the idea that conduction can set 
in only if the energy density exceeds a critical level essen¬ 
tially everywhere in the sample. However, this neglects 
the fact that local fluctuations away from the average 
energy density generally cause a breakdown of pertur¬ 
bation theory and may induce delocalization. Indeed, 
the perturbative analysis of Ref. [5[ focuses on scattering 
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processes where at every vertex an additional particle- 
hole pair is created, which was justified by the paramet¬ 
rically larger number (K) of such diagrams as compared 
to diagrams that preserve or even reduce the number 
of particle-hole excitations. However, going beyond this 
approximation and starting from a generic initial low en¬ 
ergy state, after a finite number of steps in the perturba¬ 
tive expansion one couples to configurations which con¬ 
tain hot, internally ergodic bubbles, for which perturba¬ 
tion theory is not controllable anymore. Our arguments 
following suggest that the further perturbation theory 
on this diagrammatic branch diverges and cannot be re¬ 
summed. We conjecture that a parametrically small (in 
K) but finite conductivity results in that case. 

In physical terms, we argue that delocalization occurs 
as soon as finite, but mobile excitations exist, even if 
they occur with very low density. These highly excited 
fluctuations constitute the analogs of the diffusive pairs 
in the assisted hopping model discussed above. Examples 
of such excitations are large, albeit finite regions which 
are hotter than their environment and thus are internally 
ergodic. Hereby we assume that interactions are local, 
so that the internal ergodicity is only a function of the 
energy contained in that region. 

Let us stress already now that the bubble excitations 
considered in this paper are thermal, and are not tied 
to anomalous realizations of the disorder in specific loca¬ 
tions. Therefore, those bubbles can potentially show up 
in any location in the system. We notice that, as far as 
thermal bubbles are concerned, the strategy of Imbrie [l| 
to show the existence of an MBL phase would fail. In¬ 
deed, it requires that the location of all possible resonant 
spots can be determined independently of the state of 
the system. As this is a crucial point, we review in Ap¬ 
pendix [B] the treatment of Ref. [l[ of ergodic spots that 
are tied to rare disorder realizations in specific locations. 
In particular we spell out why those do not lead to delo¬ 
calization of fully MBL systems in d = 1, while thermal 
bubbles do. 

Let us now assume that at some temperature there is 
conduction and ergodicity. [22| In typical states and in any 
given place finite bubbles of sufficiently high temperature 
(that ensure internal ergodicity) occur with finite proba¬ 
bility as spontaneous fluctuations of the energy density. 
Those are not tied to a particular local disorder realiza¬ 
tion. Thus, at any instant of time there exists a possibly 
very low, but finite density of such ergodic fluctuations. 
Following, we argue that such excitations are mobile and 
delocalize the whole system, akin to the diffusing pairs 
above. From this reasoning it follows that finite con¬ 
duction at some temperature implies finite conduction at 
any temperature in thermodynamic systems with local 
interactions. [23 As a consequence, systems in the con¬ 
tinuum should exhibit finite transport at any T > 0, as 
they always possess ergodic states at high enough energy 
(see also the discussion in Ref. 0 ). 

To argue for the mobility of bubbles, we show that 
there exists a resonant, delocalized subset of bubble con¬ 


figurations. We consider a quantum lattice system with 
local interactions and a bounded energy density, possess¬ 
ing a putative many-body mobility edge at energy den¬ 
sity Cc, such that states below (above) Cc are localized 
(ergodic). For simplicity, we assume the model to be 
one-dimensional. Now, consider a rare hot bubble of a 
super-critical energy density at some 62 > Cc, surrounded 
by ’’cold” regions of energy density ei < Cc- If this energy 
fluctuation is large enough [much larger than a correla¬ 
tion length ^( 62 )] and decoupled from its surrounding, it 
is internally ergodic by assumption. 

We argue that this state can hybridize with a translate 
of the bubble by some length io > max[^(ei), (^( 62 )] when 
the coupling between the hot region and its surrounding 
is switched on. It suffices to show that extending (or 
shortening) the hot region by a length Iq (by heating up 
or cooling down the neighboring regions of size io) can 
occur as a resonant transition. The hybridization with 
the translated bubble then follows from two successive 
hybridization processes, as illustrated in Fig. [2l For the 
latter, it is enough to show that changing the energy in 
the boundary region by a finite amount is a resonant 
process. Let Hi = gOh ^ Oc be the interaction term 
coupling a hot (h) and a cold (c) region of size ^0 across 
their common boundary. Let T, T' be eigenstates in the 
hot region and 77 , rj' eigenstates in the cold region. For 
any hot eigenstate T in a sufficiently large bubble we can 
find (many) T' such that 

|(v|/7?|gl|v|/V)| 

because on the one hand, by the eigenstate thermaliza- 
tion hypothesis (ETH) [^, |(T |0/j,|T')| ^ where 

dh is the dimension of an appropriate micro-canonical en¬ 
semble for the hot bubble at the energy density set by T, 
while the matrix element \{r]\Oc\r]')\ = 0 ( 1 ) is finite and 
independent of dh- On the other hand, we can pick T' 
such that \E{ri) — E{ri') ^ E{^) — E{^')\ < kF/d/,,, where 
W is the energy width of the ensemble. The ratio in ([ 2 ]) 
thus scales as ^ d^ and grows exponentially with the 
length of the bubble. It may thus become much larger 
than unity, indicating a resonant process. This is not 
surprising: it merely expresses that a sufficiently large 
ergodic bubble acts as a bath for small systems coupled 
to it. It follows that configurations with hot bubbles in 
different positions hybridize with each other. We expect 
that the eigenfunctions of the system hybridize essen¬ 
tially all configurations which are resonantly connected, 
implying delocalized eigenfunctions. Since it is easy to 
check that any configuration consistent with global con¬ 
servation laws can be reached via resonant processes, we 
expect that eigenstates also satisfy ETH in the thermo¬ 
dynamic limit. While for generic many-body systems our 
arguments rule out the coexistence of localized and delo¬ 
calized states, mobility edges are instead well established 
for one-particle systems. In Appendix [Cl we explain that 
this does not imply any inconsistency, since our reason¬ 
ings about bubbles do not apply to one-particle systems. 
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In Sec. imi we will present a more thorough discussion 
of the properties of the resonant subgraph of configura¬ 
tions. 



FIG. 2: Hybridization process. The state with the bubble 
on the left (top) hybridizes with the state with the bubble on 
the right (bottom) via an intermediate state. Equation (|2} 
shows that the two transitions depicted here are resonant. 


C. Potential caveats: Can bubbles freeze despite a 
percolating resonant subgraph? 

We now ask whether processes that have not been 
taken into account in the previous analysis could impede 
the hybridization of bubbles. Indeed, what was argued 
up to now can be summarized by saying that there ex¬ 
ists a connected subgraph of base states, along which all 
transitions are resonant. This subgraph is delocalized in 
the sense that most member base states differ from each 
other in degrees of freedom at arbitrary far distances in 
real space. These are the base states that contain well- 
delimited and ergodic bubbles. What happens when we 
consider the coupling to base states lying off this reso¬ 
nant subgraph? However, before we go into this techni¬ 
cal analysis [see (b) below], let us first look at it in an 
intuitive fashion. 

(a) Complete disappearance of bubbles. If one thinks 
about the issue in a dynamical setting, rather than as an 
exercise in spectral perturbation theory for eigenstates, 
then the following consideration comes up: Intuitively, 
hot bubbles should not survive dynamically, but should 
rather spread, dilute their energy, and eventually localize, 
so that they could not evolve back to their original hot 
configuration. If this were true, then surely it would rule 
out delocalized eigenstates. 

However, although such a spreading is indeed entropi- 
cally favored in real-time dynamics, the above conclusion 
that hot configurations could therefore not form again 
is fallacious. We consider an ensemble of mutually or¬ 
thogonal initial states of a given energy density. When 
summing the projectors onto these initial states, we re¬ 
cover a thermal (microcanonical) density matrix, even if 
the initial states themselves are not thermal, and hence 


the mean number of bubbles remains constant in time 
by the time invariance of the thermal density matrix. In 
simple terms: thermal fluctuations do occur, and keep 
occurring, irrespective of whether the system is localized 
or not. The only thing that localization could do is to 
pin the fluctuations at fixed positions. However, our ar¬ 
gument is ruling out precisely this option. 

(b) Technical analysis of off-resonant couplings. We 
should convince ourselves that coupling to configurations 
off the resonant subgraph does not spoil the resonance of 
transitions on the graph of bubble configurations. Sup¬ 
pose that the spreading mixes the original hot bubble 
states with states in which the bubble has spread par¬ 
tially, and which thus form a larger local Hilbert space of 
finite dimension d'^ ^ dh To argue as much as rea¬ 
sonably possible against delocalization, we suppose that 
the resonant coupling between hot regions centered in 
different positions is restricted to the original dh config¬ 
urations (because the spread-out bubbles may have lost 
their ability to translate directly). Assuming ergodicity 
within the larger space of dimension Dh = dh ^ d '^, the 
matrix elements get reduced by a factor Dh/dh (see Sec. 
lYl for calculations). However, the minimal denomina¬ 
tors decrease by essentially the same factor to WfDh^ 
albeit with a slightly larger energy range W. The ra¬ 
tio W/W is bounded, hence it cannot destroy the hy¬ 
bridization if the original bubble was large enough. This 
contrasts with single-particle problems where sufficiently 
strong coupling to a bath may induce localization due 
to a significant increase of the effective bandwidth, as 
discussed in [^ . 

In summary, the admixture of configurations that are 
not part of the resonant network cannot prevent the reso¬ 
nant hybridization along the network, but it does increase 
the timescale for transitions between different positions 
of the bubble by a factor Dh / dh . This is ultimately very 
similar as in the assisted hopping model where the possi¬ 
bility of separation of two spatially close particles cannot 
prevent their finding back together and diffusing further, 
while it can diminish the pair diffusion constant. 

(c) Other obstructions to delocalization. There are ex¬ 
amples of nonergodic behavior that are not straightfor¬ 
wardly captured by the above analysis: ( 1 ) single-particle 
localization in weak disorder in low dimension, where the 
proliferating amplitude of return to the origin in d < 2 
destroys hybridization at large distances; ( 2 ) the ’’quan¬ 
tum percolation” problem, where we see generation of 
random self-energies from the structural disorder along 
barely percolating paths [H, [29|; (3) many-body orthog¬ 
onality catastrophes, as in spin-boson systems at T = 0 
and related spin problems at finite T [^. All these ex¬ 
amples rely on specific mechanisms that seem not to be 
present in our problem. For the case (1), this is addressed 
explicitly in Sec. IIVI Another issue are rare regions with 
anomalously strong disorder. Those may render trans¬ 
port in d = 1 subdiffusive [3l|,[32|, but we argue in Sec. IIVI 
that they do not prevent delocalization by bubbles. 
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D. Analogy with bubbles in kinetically 
constrained, classical glasses 

It is interesting to draw an analogy between the bub¬ 
bles discussed here, and bubbles that mediate transport 
in classical, kinetically constrained models [33[. In both 
cases, mobility is ensured due to rare fluctuations in local 
configurations that allow the system to move and explore 
the phase space ergodically. In generic kinetically con¬ 
strained models, the diffusion constant remains finite and 
only vanishes as the system becomes entirely jammed at 
maximal density. Those models are dynamically very 
similar to disorder-free quantum models or disordered 
quantum models with a high-temperature ergodic phase, 
for which our arguments imply the impossibility of gen¬ 
uine and robust many-body localization. 

In special classical glass models, such as the Knight 
model, [ 33 I, certain moves are strictly disallowed. More¬ 
over, the possibility of moves generally depends on a 
set of configurational constraints, whose number may di¬ 
verge upon tuning a parameter. In that case, a genuine 
dynamic glass transition can take place. That transi¬ 
tion is vaguely analogous to the many-body localization 
transition, tuned as a function of disorder or interaction 
strength. Indeed, in both cases any finite bubble with 
whatever optimized properties remains ultimately immo¬ 
bile. However, in contrast to the Knight model, which 
has no quenched disorder, a genuine MBL transition does 
require quenched disorder. 

Of course, an appropriate quantum version of the 
Knight model would also be many-body localized, but the 
corresponding phase would probably not be robust with 
respect to local perturbations of the Hamiltonian since 
those re-introduce finite amplitudes for moves that were 
previously exactly suppressed. The latter constitutes an 
essential ingredient for the sharp glass transition in the 
kinetically constrained models. Insofar one may consider 
the kinetically constrained glass phases as fine tuned (re¬ 
lying on certain exactly vanishing transition amplitudes). 
On the other hand, the robustness of MBL with respect 
to random local perturbations of the Hamiltonian has no 
classical analogy to our knowledge. 


III. DETAILED ARGUMENT FOR 
HYBRIDIZATION BY BUBBLES 

We describe here our argument for hybridization bub¬ 
ble configurations in eigenstates at a more formal and 
detailed level. We consider a quantum lattice model with 
local interactions, having a putative many-body mobil¬ 
ity edge. For concreteness, we assume that the model 
is one dimensional and that states below energy density 
Cc are (putatively) localized, whereas those above Cc are 
ergodic. We choose the energy density of the bottom of 
the spectrum as a reference and set it to zero. We also 
impose a maximal energy density > Cc, reflecting the 
fact that the Hilbert space is locally finite. 


For our argument it is important to have states at dis¬ 
posal that are clearly ergodic or localized in a given finite 
volume. Therefore, we introduce, somewhat arbitrarily, 
a window of energy densities Cd < Cc < 6 c 2 that are too 
close to critical to be identified clearly as localized or de¬ 
localized in nature. Let ^(e) be the localization length, 
which diverges as e ^ Cc- We express lengths in units 
of lattice spacings, and energy densities e as energy per 
site. We will now coarse grain the model and group ^0 
adjacent sites into ’’grains”. The length is chosen such 
that it is (a) larger than the localization length ^(cd), 
and (b) large enough so that the interaction energy be¬ 
tween two neighboring grains is small compared to ^o^d, 
the maximal energy in a localized grain. The second con¬ 
straint ensures that the interaction of a low-energy grain 
(e < Cd) with its surroundings does not trivially suf¬ 
fice to render the grain ergodic, and, similarly, that the 
interaction of a high-energy grain with its surroundings 
does not trivially suffice to localize that grain. This will 
be satisfied by the choice of a large enough £ 0 , since the 
interactions are local. 

The coarse graining provides a useful starting point, 
from which to proceed with perturbation theory in the 
inter-grain coupling. Within each grain, we compute the 
eigenstates which come in three kinds: cold (below ecj, 
hot (above 602 ), or intermediate (between Cd and 602 )- 
If we consider the Hamiltonian without the interaction 
between grains, then obviously the eigenstates are prod¬ 
ucts of grain eigenstates. Let us focus on eigenstates at 
very low energy density e <C Cd- Our aim will be to show 
that even those low-energy states are delocalized. In such 
states, typically, non-cold grains appear only with a den¬ 
sity u that tends to 0 as e/cd ^ 0 . Chains of labeled 
grains such as ccciccchhhiiccc serve as ” mesostates”, 
with c/i/h standing for cold/inter mediate/hot. We can 
now write our model as 

H = Y,[Ho{x)+H^{x,x + 1)], (3) 

X 

where x labels the grains, Ho{x) acts on the Hilbert space 
1-L{x) at grain x only, and Hi describes the coupling be¬ 
tween neighboring grains. Mesostates are collections of 
eigenstates of the term Hq. We now consider switching 
on the coupling terms and evaluate their effect on the 
unperturbed eigenstates of Hq. This procedure is similar 
to the one followed in Ref. Q. First we add the inter¬ 
action terms between cold grains (see following for what 
is meant precisely). Since we assumed that ^(cd) < £ 0 , 
this will not have much effect on the localized eigenstates, 
which thus remain close to products. Note that by doing 
this, from the point of view of a typical state at energy 
density e, we have already added most of the interaction 
terms. What remains is a small fraction ^ 2u (which is 
controlled by the overall energy density) of all interaction 
terms. We now add the interaction terms between hot 
grains. By assumption, sufficiently long stretches of such 
grains ... hhhh ... (which we call ’bubbles’) are ergodic 
and we will assume that the resulting hot eigenstates in 
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those bubbles satisfy the eigenstate thermalization hy¬ 
pothesis (ETH). The situation at this moment is hence 
that we have partitioned the Hilbert space into a big di¬ 
rect sum, and the Hamiltonian Hq is block diagonal, with 
the blocks labeled by mesostates. Let be the projec¬ 
tor that restricts the value of Hq {x) so that grain x is of 
type r = h^i^c. With this notation the interaction terms 
that have already been added are 

p:- plix Hi{x,xp i ) p :^ piir. ( 4 ) 

for {rx,rx+i) = = (c,c), and for (ra,,r^+i) = 

= (h, h). Some terms are obviously very small 
(because the interaction is local in energy) and seem 
irrelevant, namely those corresponding to, {rx^Vx+i) = 
(c, c), = (h, h) and with primes and no primes 

reversed. The main terms that we will be focusing on are 
those that allow bubbles to spread and move. Those are 
terms with 

Tx = r'x = and arbitrary (5) 


a state ^ corresponding to the mesostate ccchhhhccc and 
such that Ho{y) is not far below ed- We will argue that 
^ hybridizes with a lot of states corresponding to the 
mesostate ccchhhhicc where r'{y) = i. If instead Ho{y) 
is far below ed, then it hybridizes with a lot of states 
T' which still corresponds to r'{y) = c {ccchhhhccc)^ but 
now with Ho{y) a bit closer to ed- (A direct step to an 
intermediate state might instead require adding too much 
energy in one transition. Therefore, we split the process 
into several small heating steps to make sure our argu¬ 
ment remains valid also if the interactions are assumed to 
be strictly local in energy.) Finally, we need to increase 
the energy stepwise from i to h at grain y. The argument 
for all these transitions is essentially the same and for the 
sake of simplicity, we stick to r{y) = r'{y) = c. The next 
subsection shows the flexibility of the argument. 

Obviously, it suffices to take eigenstates in in the 
region [x^y] because of the essential product structure 
(exact at the left edge, approximate at the right edge 
around because we have already included the coupling 
between cold regions). They are of the form 


and with x and x 1 reversed. They are at the focus of 
the next section. 


A. Resonant delocalization of bubbles 

Let us now consider states of the following form (hot 
bubble in a cold environment) 

cccccccc hh . . . hh cccccccc, (6) 

n grains 

where n is sufficiently large so that the eigenstates in 
the bubble satisfy ETH. We now argue that this state 
hybridizes with translates of the bubble when we add 
some of the missing coupling terms: In particular, we 
want to admix the mesostates (with x, y labeling specific 
grains) 

... ccchhhhccc ... ... cccchhhhcc ... (7) 

X y X y 

in which the bubble has been translated by one grain. 
More precisely, we mean that most microstates (i.e., 
eigenstates of the Hamiltonian considered up to now) 
corresponding to the left mesostate can hybridize with a 
lot of microstates corresponding to the right mesostate. 
This in turn strongly suggests that we should expect 
all eigenstates to delocalize completely over these two 
mesostates. To obtain this, we have included the relevant 
coupling terms (|4)) corresponding to two bonds (x, x -h 1) 
and {y — l^y). This hybridization process can be broken 
down into elementary steps, that is, transitions at first 
order of perturbation theory. First, by energy exchange 
with the hot region, the cold (c) grain at y is heated until 
it becomes intermediate {i) and finally hot {h). Second, 
the h grain at x is cooled down until it becomes c, via 
intermediate stages of i. Microscopically, let us consider 


^ = ^(g)7/, ^' = T'(g)V, (8) 

where 77 , 77 ' are the unperturbed eigenstates at grain y, 
while T, T' are hot bubble states in the region [x, y — 1] 
consisting oi n = y — x grains. Consider T' such that its 
energy (evaluated with Hq) is within a range W ^ em 
of the energy of T. The space spanned by such states 
has dimension dh ~ exp[ 5 ^ 0 ^] which grows exponentially 
in 71, 5 being the corresponding entropy density. Write 
Hi{y — l^y) = gOh 0 Oc^ the first factor acting on ^ — 1, 
the second on y. Assuming ETH, the off-diagonal matrix 
elements of local operators are given by 

(9) 


In other words, the (non-eigenstate) vector is essen¬ 
tially a random amplitude superposition of eigenstates 
T'. Take now /S.E := E{if]) — E{r]') sufficiently small, i.e. 
not exceeding W, then \{r]\Oc\rj')\ ^ 1. In fact, assur¬ 
ing the non-vanishing of | (77|Oc|^01 main reason to 

choose AE sufficiently small. We can then find many T' 
(in fact, ^ y/d^ of them) such that 


|(T77|i7i|TV)| 

|AE;-fE;(T) -e;(t')| 


( 10 ) 


because the energy spacings are of order W/dh and 
(T77|iLi|T'77') ^ Hence the ratio in m is huge 

since dh grows exponentially in n. 

The outcome of the above calculation should not come 
as a surprise: it merely expresses that an ergodic bub¬ 
ble can act as a bath for a small system (here grain y) 
that is coupled to it. Upon repeating the same calcula¬ 
tion a few times, one easily convinces oneself that states 
with the bubble in different positions hybridize with each 
other. This in turn implies that they should appear with 
comparable amplitudes in typical eigenstates. 
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B. Spatial range of direct hybridizations 


In the above derivation, we focused on hybridizations 
that result in the translation of a bubble by one grain. 
One might worry that this is too negligible a transla¬ 
tion if the bubble is very large, n ^ 1. However, here 
we show that direct hybridizations can take place at dis¬ 
tances which are a finite fraction of the bubble length. 

As already pointed out, in the above derivation, we 
were careful to pick states 77 , 77 ' whose energy difference 
was small enough so that 1(77, OcV) I ^ 1 - This is, how¬ 
ever, not crucial, and if r'{y) = then it cannot be 
assured anyhow. The matrix element | (77I Oc I 
typically decay exponentially in the energy difference 
^( 77 ) — Hence, it can be as small as but 

obviously this number decreases with £q but not with 72 , 
so it cannot compete with the latter, if the bubble is 
sufficiently large. To determine at what distance direct 
hybridizations are possible, we proceed as follows. In¬ 
stead of making the transition 77 ^ 77 ' at grain 77 , we now 
make a transition 77 ^ 77 ' in a stretch of £ grains starting 
at y. By the structure of localized states, we know that 

\{R\Oc\rf}\^ {g/em)“°- ( 11 ) 

The transitio n is po ssible as long as this small number is 
larger than so that we find 


£r. 


s 

-77/. 

2\n{em/g) 


( 12 ) 


Sec. IIV A 11 we argue that a scenario in which bubbles 
disappear dynamically is inconsistent. This is a concep¬ 
tual point. Then, in Sec. IIV A 2 [ we examine the hy¬ 
bridization argument on a formal level and we exclude 
that bubbles get localized by quantum dynamical effects. 
In Sec. IIVBI we discuss how the fact that bubbles tend 
to spread slows down their motion. 


1. Persistence of bubbles 

Let us write (•)£ for the expectation value in a mi- 
crocanonical ensemble at energy density e (containing a 
large but subextensive number of eigenstates). Let n be 
the minimal length (in units of grains) of a well-ergodic 
bubble, as considered above. The thermodynamic prob¬ 
ability of having such a bubble around position x is given 

by 

Px — {Ex)ei Ex = X{E[x—n/2,x^n/2] ^ ^c2'^£o)i ( 13 ) 

i.e. Bx is the indicator of a local fluctuation at sufficiently 
high energy density. Of course, Px becomes very small as 
ec 2 — e, or 77 /, or both are taken large, but it remains finite 
and independent of the total volume. By definition of 
the microcanonical ensemble, 

= ( 14 ) 


This shows that the bubble hybridizes directly with a 
bubble configuration translated by a finite fraction of its 
size. However, this fraction becomes parametrically small 
as the coupling becomes weak ^ 0 . 

IV. DISCUSSION OF POTENTIAL CAVEATS 

We now discuss in more detail the caveats introduced 
in Section iHCl 

A. Can bubbles spread and permanently localize? 

We have not yet added all coupling terms from Hi. 
Indeed, not only can the bubble move through the cold 
background, it can also spread its energy. Entropically, 
this is of course much more likely in real-time dynamics. 
In particular, one sees that starting from a bubble con¬ 
figuration, the most likely thing to happen dynamically 
is that the bubble spreads until its energy density is in¬ 
termediate or just below the putative mobility edge. At 
that point we cannot expect it to spread further as the 
involved states are now localized. The question arises as 
to whether these further couplings may induce a local¬ 
ization of bubbles, despite the above construction of an 
apparently resonant, delocalized network of bubble con¬ 
figurations. We address this issue in two steps. First, in 


and this means that it is not possible that no eigenstate 
has any weight on bubble configurations. One possibil¬ 
ity consistent with da is that typical eigenstates have 
an appreciable weight, on bubble configurations at 
X. Another (extreme) possibility is that a fraction px 
of eigenstates have weight nearly 1 on bubble configura¬ 
tions, while the others have none. Either way, this rules 
out the scenario that bubbles could be completely absent 
in the eigenstates. 

One can also formulate the persistence of bubble con¬ 
figurations in a dynamical way. This is more appealing 
if one thinks of our arguments as pertaining also to the 
evolution of well-chosen initial states, that contain a def¬ 
inite bubble. Let us show this in a special and simple 
case. Consider a model where particle density is a con¬ 
served quantity and there is a putative mobility edge as a 
function of density, even at infinite temperature. Two ob¬ 
vious examples are the models studied in Refs. mu. In 
such a case, a relevant ensemble is one which constrains 
the particle density ^ G [p — 6^ p -£ S\ to be close to p, 
without any constraint on energy. The advantage of this 
case is that the projector onto the equilibrium ensemble 
corresponding to this constraint, can be decomposed into 
a basis of initial states |s) which are products over grains 
with a definite particle number on each grain, such that 
the total density is indeed in [p — 5^ p P 5]. Now, this en¬ 
semble |<^)(«5| is exactly invariant under the dynamics 
and hence we obtain that the expectation value for seeing 
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a bubble at position x, 

'^{s\UlBxUt\s) = Px, for any t, (15) 

is invariant under time evolution, and remains finite at 
all times. Here and Px are again the indicator of 
a bubble around x and the thermodynamic probability 
of a bubble (defined analogously to (p!3|) ). and Ut is the 
Hamiltonian time-evolution. 


2. Robustness of hybridization 


In the notation of section HTTl the relevant type of tran¬ 
sitions are the following: 

ccchhhhccc ^ cccchhhhcc 

X X 

... ... cciiiiiiiicc ccciiiiiiic ... ..., 

(16) 

where the states on the lower line represent a multitude 
of mesostates in which the bubble has partially spread its 
high energy density. Let us assume that those states do 
not communicate with each other. This simplifying as¬ 
sumption favors maximally the possibility that the cou¬ 
pling to such states could localize the bubble and thus 
could invalidate our preliminary conclusion regarding de- 
localization. We now consider the two subspaces, each 
of dimension that correspond to the mesostates on 
the upper line, the eigenstates of which are hybridized 
by the perturbation Hi. Let us refer to them as left and 
right subspaces. We then couple each of them to a space 
of spread bubbles, having a dimension d'^ ^ dy^ and ask 
whether the perturbation Hi is still able to induce hy¬ 
bridization between left and right subspaces. Concretely, 
the subspace is now embedded in the space 
of dimension Dy = d'^-\- dy, and the inter-grain coupling 
operator Oy becomes 0/^00. We focus on the transitions 
between the ergodic states T, T' (notation as above), and 
just consider the operator Oy which acts on the hot bub¬ 
ble. Let us assume that after diagonalizing within the 
larger spaces of dimension the eigenstates T, T' are 
completely ergodic within those spaces and well captured 
by random matrix theory. (In practice, this defines the 
relevant space to which the bubble subspace should be 
extended, and its dimension Dy.) We now have to dis¬ 
cuss how the ratio 


\AE + 

differs from the original ratio 


(17) 


\AE + E{^) - Ei'^')] W ^ ’ 

with given |AL^| < W. We find a suppression of the 
numerator because now 




Dh ■ 


(19) 


Indeed, the simplest way to derive this is by remarking 
that 


^ |(^|0^|^')|2=IY(0j.0^)~4, (20) 

as Oh acts only in the original subspace (with dimen¬ 
sion dh) and it is zero on the attached space with di¬ 
mension d'^. On the other hand, the energy spacing 
\AE-\-E{^)—E{^')\ can now be made as small as W/Dh, 
where W is the width in energy of all states that signif¬ 
icantly couple to the original bubble states. It follows 
that the ratio (HZI, and hence m, is reduced by a fac¬ 
tor W/W. If this effect rendered the ratio (pT|) smaller 
than 1, the eigenstates would likely not hybridize across 
the subspaces, i.e., we would find localization induced by 
coupling to further degrees of freedom. However, the max¬ 
imal conceivable value of W is of order with iy the 

length of the region to which the energy spreads. En¬ 
ergy conservation and localization below Cc lead to the 
upper bound ih{^c — e) < n{em — e) (recall that e < Cc 
is the typical energy density in our system). This yields 
W/W ^ C with C independent of n. This is insufficient 
for localization, since the ratio m is exponentially large 
in n. Thus, the hybridization of states with large enough 
bubbles survives, despite their spreading to entropically 
more favorable states. This contrasts with single particle 
problems where the coupling to extra degrees of freedom 
was found to induce localization under certain circum¬ 
stances BUI. In those cases, there is no exponen¬ 
tially large factor that offsets the effect of an increased 
bandwidth W, which renders coupling-induced localiza¬ 
tion possible. 


B. Dynamic retardation 

Even though the inclusion of the states on the lower 
line of (pT|) cannot prevent hybridization, it does of course 
increase the timescale necessary for transitions between 
the two bubble positions. The transition rates can be 
estimated from a simple Eermi Golden Rule calculation 
as 


1 

bef \AE + E{^) - E{^')y ^ ’ 

|AE + E(5'')-^(5')!’ 

before and after including the extra states, respectively. 
The first rate is of order g^/W, while the second is of 
order [dy/Dh)g^/W. Hence, by adding the new states, 
we have increased the timescale by order Dy/dh (keep¬ 
ing only terms exponential in n). This is very intuitive: 
Transitions are now only possible from a fraction dh/Dy 
of all states, and accordingly it takes longer until a transi¬ 
tion will be attempted. Alternatively, one can view this 
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as follows: For a large bubble close to criticality (with 
structure cciiiiiiiiicc) the ’active’ configurations of the 
type cccchhhhcccc manifest themselves as large devia¬ 
tions, which occur with exponential rarity (in the bubble 
size). Yet, as shown above, they do lead to percolating 
hybridization of eigenstates, and hence to delocalization. 

C. Harmlessness of rare, 
strongly disordered regions in 1 dimension 

So far we have tacitly assumed that the existence of a 
global thermal phase at some energy density implies that 
any large enough finite region has ergodic states at that 
energy density. However, due to rare fluctuations of the 
disorder, it can happen that rare, large regions still have 
all their states localized as long as they are disconnected 
from the rest of the system. One may think that espe¬ 
cially in d = 1 such regions could block global transport 
and thus localize the system at low temperature. Even 
though this effect further increases the time scale neces¬ 
sary for thermalization, we now argue that it does not 
prevent it. 

As a preliminary, we consider a fully localized system 
of length Lq in contact with an ergodic system of length 
Li . By an analogous argument as used above for resonant 
delocalization, we see that for LofLi smaller than some 
number [depending on the localization length and the 
entropy density of the ergodic system, see e.g. Eq. m] 
the coupled system will be ergodic: all formerly localized 
states can hybridize with each other. 

Now to the main argument. Let us consider ergodic 
bubbles of some large size £. Now consider rare regions 
of exceptional disorder that could block such bubbles, 
since an adjacent bubble of size i would not be able to 
heat up this region to ergodic states (of type h). Let 6 = 
S{£) be the typical distance between such rare regions. 
By the above preliminary remark, such blocking regions 
have a length ^ £. However, being rare regions (large 
deviations), the typical distance between them is much 
larger, 6{£) for some c > 0. Now consider a bubble 

between two blocking regions. It renders the whole region 
between them ergodic. Hence, the blocking regions are 
in fact next to an ergodic bath of length S{£), which is 
exponentially large in £. Accordingly, its level spacing is 
double-exponentially small in £. Thus, tunneling under 
the barrier of thickness ^ which is only exponentially 
small in £, will easily hybridize the ergodic regions on 
either side and ensure transport. 

D. Bubbles and weak localization in low 
dimensions 

As is well known, in dimension d < 2, non interacting 
particles are weakly localized by any non-vanishing disor¬ 


der strength. Since bubbles resemble a particle-like exci¬ 
tation, we should discuss whether they undergo a similar 
weak localization in low dimensions. 

Indeed, at zero temperature, the answer is expected 
to be positive. However, since we consider a finite en¬ 
ergy density, it would be incorrect to picture the bubble 
as moving in a low-dimensional fixed disorder potential. 
As the bubble moves, it can excite or relax degrees of 
freedom. Thus, the Hilbert space locally resembles a 
tree, rather than a low-dimensional lattice (the number 
of relevant configurations that can be reached as the bub¬ 
ble moves grows exponentially, rather than polynomially 
with the traveled distance), and thus weak-localization 
effects should become irrelevant. 

More concretely, let £ be the localization length of the 
bubble motion at zero temperature, i.e., in the ground 
state (fixed environment). Obviously, £ increases with the 
bubble size since larger bubbles have more internal states 
(this is analogous to the increase of the single particle 
localization length with the number of channels in d = 
1); for large bubbles £ will be due to weak-localization 
effects. Now, consider finite energy density, and d = 1 
for simplicity. Let 5 > 0 be the entropy density in the 
cold background. Then, the condition 5-^ 1 is sufficient 

to ensure that inelastic scatterings of the bubble occur 
before the weak localization manifests itself, rendering 
them irrelevant for large enough bubbles. Note, however, 
that this condition places an additional lower bound on 
the size of mobile bubbles in low dimensions. 


E. Restrictions on the type of rare events that may 
lead to delocalization 

The reader may wonder whether certain rare events 
similar to the bubbles discussed here would not rule out 
the possibility of genuine many-body localization alto¬ 
gether, or (erroneously) imply the absence of mobility 
edges in single-particle cases as well. However, in two 
appendixes we explain that this is not so. 

In Appendix [B1 we contrast the bubbles discussed here 
(namely disorder-independent high-energy fluctuations) 
with rare spots of low disorder in an otherwise fully lo¬ 
calized system. We show why the latter are benign and 
do not hamper Imbrie’s approach [l[ to demonstrate lo¬ 
calization, whereas high energy fluctuations do destroy 
localization, if the system is ergodic at high tempera¬ 
ture. In Appendix [Cl we explain why our considerations 
apply only to genuine many-body systems, and how the 
bubble construction fails when applied to a single-particle 
system. 


V. NUMERICAL RESULTS 
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x/L x/L 

FIG. 3: Left: Disorder averaged energy per link Si ai t = 0 (red) and averaged over time (green) for L = 12. Initially a cold 
region of length Lc = L/2 is prepared. The disorder strength is SJ = 3J. Right: Same protocol, but for SJ = J and very 
short cold intervals {Lc = 2), at various L. The memory effects diminish with increasing L, but the hot region fails to 
thermalize the system well, even at the largest sizes. Results were averaged over 5000 disorder realizations. 


Our theoretical arguments contradict recent numeri¬ 
cal data in favor of mobility edges The inconsis¬ 

tency is, however, only apparent. Indeed, we find that 
numerically accessible system sizes are not sufficiently 
large to host bubbles that are ergodic enough to be mo- 

I 


We study the disordered Ising chain with next-to- 
nearest neighbor interaction considered in Ref. Q, 

L 

H = - + 6Ji) crK+i + J2<0-f+2 + , 

i=l 

where SJi G are independent random vari¬ 

ables, and periodic boundary conditions are taken. We 
choose parameters J = 1, J 2 = 0.3 and = 0.6 as in 
Ref. 0, but add a finite = 0.1 to remove the Ising 
symmetry and the associated degeneracies. The phase 
diagram in Ref. Q predicts a mobility edge in the ther¬ 
modynamic limit at disorder strength SJ = 3. To test 
our ideas, we prepare the system at = 3 in a product 
state ofthe form |'0(O))^ = {(fc)L,^\Xh) l-l^^ l^c) 

is the ground state of an interval of Lc sites, while \xh) 
is an eigenstate of the complement close to the middle of 
the spectrum (a hot bubble). We choose L—Lc as large as 
possible but such that the resulting global energy density 
is below the putative mobility edge. We then compute 
the time-evolving energy density on link (i, i -h 1): 

Si (t) = - (J + SJi) (?/) (t)| crfcr|+i (t)) . (23) 

Our theory of mobile bubbles would predict that the Ci {t) 
profile becomes approximately fiat as t ^ 00 . Via exact 
diagonalization, we evaluated its time average in finite 
system sizes, but almost no energy spreading from the 


bile. Therefore, delocalization by bubbles could not have 
been seen in numerics up to now. In other words, the nu¬ 
merical results do not contradict delocalization by rare 
bubbles, but rather confirm that available sizes are not 
large enough. 


initial state was observed [cf. Fig. [3] (left)]. For tiny cold 
regions (Lc = 2) and hot ’’bubbles” of almost the system 
size the global energy density is supercritical. Yet, still 
only a very small fraction of the bubble energy spreads 
to the cold region at L = 12 (not shown), while in the 
thermodynamic limit, the energy profile would obviously 
thermalize and become fiat, as the system would be in 
its ergodic regime. Therefore, these data show unam¬ 
biguously that at our system size the hot region is still 
unable to act as a bath. 

To document this further, we considered normalized 
states af ja), with a an eigenstate, and calculated the 
inverse participation ratio (IPR) of its decomposition into 
eigenstates |/3) of the full system: 

IPR„ = ^|(^|aJ|a)G (24) 

h 

The results are shown in Fig. El At strong disorder, 
eigenstates are nearly eigenstates of erf as well, and thus 
IPRck ~ 0(1), with a distribution expected to become 
system-size independent for large L. As shown in the 
bottom-right panel of Fig. El this is indeed the case for 
very strong disorder SJ = 5J, for which Ref. 0 found 
the whole spectrum to be localized. Conversely, deep 
in the delocalized phase, one expects eigenstate ther- 
malization and behavior akin to random matrix the¬ 
ory, I (/31 erf \a)\ (X exp[—sL/2], leading to a typical value 
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6J/J=0.1 6J/J=1 



6J/J=3 


6J/J=5 



FIG. 4: Distribution of In(IPR) associated with matrix elements of al evaluated on eigenstates randomly picked from the 
middle of the spectrum, for 5J/J — 0.1,1, 3, 5. In the ergodic phase, the typical IPR is exponentially small in the size L. In 
the localized phase, the distribution is size independent. At SJ = J and the considered L, the typical IPR is exponentially 
small as in a truly ergodic phase, but the distribution is broad, and shows a tail towards localized values. SJ = 3J is nearly 
critical: the ’localized’ peak at IPR = 0(1) slowly decreases with increasing L. 


IPRc, ^ exp[—sL], with a narrow distribution. For very 
weak disorder, 6J = O.IJ, we found the exponent to be 
So. 1 ~ 0.55 which equals essentially the thermal entropy. 
The standard deviation of the distribution scales with 
system size in the same fashion, but is approximately 10 
times smaller than the mean value. This is seen in the 
data of the left-top panel of Fig. 01 where the relatively 
sharp peak in the distribution of log(IPR) has an essen¬ 
tially I/-independent width. 

Let us now discuss close to critical disorder, SJ = 3J: 
the results shown in Fig. 0] confirm the absence of a truly 
ergodic phase up to L = 12. In fact, the distribution 
of IPR’s at these parameters looks more characteristic of 
localization. Nevertheless, a slight, but clear tendency 
towards enhanced delocalization with increasing size is 
seen. This hints that in the thermodynamic limit the 
system will become ergodi^ in agreement with the finite 
size extrapolation in Ref. Q • To chart the lack of ergod- 
icity at small sizes, we also look at 6J = 1, where Ref. Q 
suggests that most eigenstates are delocalized, even at 


L = 12. Nevertheless, here, too, we find strong deviations 
from fully ergodic behavior, using the same two protocols 
as above. Even in the extreme case of Lc = 2 in Fig. [3] 
(right), despite some energy transfer, the hot and cold 
regions are still clearly distinguishable after a long time 
evolution. To quantify this effect, we consider the time 
average of the energy imbalance between hot and cold re¬ 
gions, Ae = where c denotes 

the single link fully in the cold region. The imbalance 
decays exponentially with system size, ^ exp (—L/^) 
where ^ increases with disorder strength. For 6J/J in the 
range [1,1.5] we estimate ^ ~ 0(10) (see below), which 
sets a characteristic scale required to observe genuine er¬ 
godic behavior. This suggests strongly that at reachable 
sizes the hot bubble is far from being ergodic. 

Also Fig.0]illustrates that systems with 6J = J^L = 12 
are far from the thermodynamic limit: the distribution of 
In (IPRc^) is much wider (as compared to the mean) than 
in a clearly ergodic sample. To be more precise, the mean 
value still scales exponentially as IPRc, ^ exp[—siL], 
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with Si ^ 0.5, but the standard deviation is much larger 
than in the ergodic phase: indeed, the ratio between the 
standard deviation and the typical value is of order 0(1). 
This means that there is a finite probability for a bubble 
to find itself in a position where the coupling to neighbor¬ 
ing regions does not allow it to act as a bath. A consid¬ 
erable fraction of eigenstates will thus be localized. This 
illustrates that ETH, and hence our bubble arguments, 
cannot be applied for = J, L = 12. 

Finally, we studied the time average of the energy im¬ 
balance as a function of system size and disorder 
strength 6J. We considered an initial state with a region 
of length Lc = 2 in its ground state and the remaining 
system in an eigenstate near the middle of the spectrum. 
The results are shown in Fig. [5l where each plot corre¬ 
sponds to a different disorder value. The results were 
averaged over 5000 disorder realizations. The data show 
that the imbalance decreases only slowly with system 
size. The dependence is consistent with an exponential 
decay. We have fitted the associated characteristic length 
^< 5 j, which grows with increasing SJ. These lengths ^sj 
are of the same order as system sizes achievable in current 
numerical studies. Therefore, we conjecture that in order 
to be able to observe bubbles acting as good baths, one 
would need to study systems, which are larger by several 
^< 5 j’s. Perhaps the density matrix renormalizati on g roup 
(DMRG) techniques used recently in Refs. [3J, can 
shed new light on this. 

In summary, the numerical analysis provided in this 
section clearly shows that the available system sizes are 
too small for ETH to be safely applied. Therefore, they 
are outside the range of applicability of our bubble ar¬ 
gument. The numerical data present in the literature at 
the moment thus do not disprove our considerations. 

VI. CONCLUSION 

We have argued that in the thermodynamic limit 
many-body-localized and ergodic states cannot coexist, 
not even at very different energies. This has important 
consequences on the nature of the MBL transition. On 
a lattice, it implies that a transition is possible at best 
upon tuning the interaction strength, but not the tem¬ 
perature. In the continuum, genuine MBL is replaced 
by a strong crossover in the conductivity instead, which 
is notoriously hard to distinguish from a genuine transi¬ 
tion. Nevertheless, for practical purposes, on very long 
time scales, such badly conducting phases will behave as 
if they were genuinely many-body localized. 
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Appendix A: Assisted hopping model 

Here we describe our numerical analysis for an assisted 
hopping model. The main aim is to show that delocal¬ 
ization on a resonant subgraph remains robust to adding 
additional terms that connect that subgraph to localized 
states. We also show coexistence of localized and de¬ 
localized states, a failure of Mott’s argument, which is, 
however, a particularity of the zero density limit of the 
considered model. 

Description of the model To reach the largest pos¬ 
sible system sizes, we consider a Hamiltonian in d = 2 
with spin-orbit coupling, which gives rise to weak anti¬ 
localization and thus allows for a genuine delocalized 
phase. To the best of our knowledge, this is the smallest 
system where delocalization can be expected, and is thus 
best suited for a numerical analysis. Here, ’’smallest” 
means that the dimension of the Hilbert space grows at 
the slowest possible rate with growing linear size L. 

Let H be the Hamiltonian of two indistinguishable 
hard-core bosons (with positions qi^ 2 ) having a single 
spin-^ degree of freedom, s, attached to them. We con¬ 
sider points q = (x^y) on the lattice (fLjV)^ and we im¬ 
pose periodic boundary conditions. The full Hamiltonian 
is 

= + (Al) 

where ido is the uniformly distributed on-site potential 

i^o = F W. (A2) 

Q 

Hi is the single-particle hopping Hamiltonian 

Hi = H-UgUg/), (A3) 

q^q' 

{q ^ q' denoting nearest neighbors) and H 2 is the assisted 
hopping, including a spin-orbit interaction. We describe 
H 2 by its matrix elements. Let 

*5 = {gi = (a;i,yi), 52 = (2:2,^2) : qi + 92, 

max{|a:;i - a:; 2 |, I 2/1 - 2 / 2 I} < 1} (A4) 

be the set of pairs of spatially neighboring points. We 
then define (g^, ^ 2 , «§) to be 

^siq'i,q2)^s{qi,q2) {q[,q2,s'\Hso\qi,q2,s), (A5) 

where the characteristic functions Is ensure that the ini¬ 
tial and final pair configuration belong to S. Further, 
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FIG. 5: Energy imbalance Ae as a function of system size, for disorder strengths J < SJ < 1.5J. The curves are simple hts to 

exponentials. 


where (•) is the disorder average and AE^ = ^n+i ~ 
with En being the ordered eigenenergies of the system. 
(A6) For the three classical ensembles, they take the values 


Hso = with 


= 


(x)rp _ ^{y)rp 

yi ^ X 


(cr(®) - (j(y)) 


T T — 

-‘-Xi-Lyi 


ri^) 


T^y)) 


T 

y-Xi-Ly^ 


H.c. 


Here, are Pauli matrices acting on the spin degrees 

of freedom, while the translation operators are defined 
T:ci\{xi,yi),{x 2 ,y 2 ), s) = \{xi + l,yi),{x 2 ,y 2 ),s) and 
similarly for Ty ^. is defined analogously. 

The Hamiltonian H^q is a lattice version of the Rashba 
Hamiltonian ~ We notice that restrict¬ 

ing the definition of Ffgo the first term —— 
a^y^Tx^} would lead to a degeneracy due to the lattice 
structure. This would prevent E[ from being a generic 
GSE Hamiltonian for any value of /i 2 - 

Numerical results. In all the simulations, we take L = 
9 and W = 1. The analysis is divided into two parts: 

(i) Delocalization via assisted hopping. First we take 
hi = 0 and h 2 > 0 (only assisted hopping). Since the 
majority of states (all configurations outside S) are now 
trivially localized, we restrict ourselves to the subspace 
Es spanned by all the classical states in S (see dSi), 
each coming with spin up/down. We aim at finding h 2 
such that Hq + /i 2^2 can be considered a “typical” GSE 
matrix with truly delocalized eigenstates. For this, we 
evaluate numerically the parameter r defined as 

dim{l-Ls) - 2 ^ maxjAEn, AEn-i} / ’ 

(A7) 


r(GOE) - 0.53, r(GUE) - 0.60, r(GSE) - 0.67. 

(AS) 

For h 2 = 0.7, we find r = 0.64 ± 0.05. This value is 
significantly larger than r(GUE). The discrepancy with 
r(GSE) presumably arises from the contributions from 
the more localized edges of the spectrum. 

To characterize (de)localization we use the logarithm 
of the inverse participation ratio, 

loglPR(tA) = - logio ( El I 1^) ’ 

r] 

where the sum over 77 runs over the classical particle con¬ 
figurations. 

Note that dim(H 5 ) = 648, and thus logIPR('0) ^ 2.5 
for a fully delocalized state 7 /^. From the point of view of 
the parameter r, /12 = 0.7 is rather optimal: The spec¬ 
trum is mostly delocalized, but the Hamiltonian is still 
genuinely GSE. Indeed, when h 2 becomes significantly 
larger than 0.7, the localized tails of the spectrum are fur¬ 
ther suppressed, but the value of r starts bending down 
as an effect of approaching the integrable limit h 2 ^ 00 . 

(ii) Robustness of delocalization against addition of 
single particle hopping. Let us now fix h 2 = 0.7, but vary 
hi > d. We determine numerically the statistics for the 
logIPR’s of the eigenstates 7 /; of H. The results are shown 
in Fig. [6l The central message of that data is the fol¬ 
lowing: Adding a finite hi, which connects the resonant 
subspace S to its much larger localized complement, does 
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FIG. 6: Statistics for the logIPR’s of all eigenstates of H for /i 2 = 0.7 and different values of hi. From left to right: 
hi = 0.01, hi = 0.07, hi = 0.15. Averages are taken over 500 realizations. 


not destroy the delocalization on the resonant subspace, 
as shown by the left and middle panels of Fig. [6l In par¬ 
ticular, for hi = 0.07 (middle) we see delocalized states 
(inside the subspace Hs) coexisting with a majority of 
localized states. Obviously a relatively large hi leads 
to delocalization of almost all states, with logIPR’s that 

I 


Appendix B: Robustness of full MBL against the 
presence of localized ergodic spots 

It is instructive to compare the situation of high energy 
bubbles discussed in the main text with rare spots of 
low disorder in an otherwise fully localized system. In 
particular we explain here why the latter do not induce 
delocalization, even though a naive adaptation of some 
arguments we used in the main text might suggest so. 

For simplicity we discuss systems in d = 1. In the 
above reasoning we used the fact that, if an ergodic sys¬ 
tem of length I/i is put in contact with a small localized 
system of length Lq and if Lq/Li is smaller than a certain 
number, then localization does not persist: all formerly 
localized states hybridize with each other by mediation 
of the ergodic system they couple to. Now, one could 
naively think that this argument could be iterated ad 
infinitum: if one assumes that all states in the joint re¬ 
gion of length I/O + 1/1 are now fully ergodic, so that ETH 
holds, one could try to iterate the previous argument and 
conclude that yet a larger region would become ergodic, 
and so on. This would lead to the wrong conclusion that 
MBL systems simply cannot exist, because there is al¬ 
ways a finite probability of ergodic inclusions of some 
finite size Li. However, a more careful analysis reveals 
where this reasoning is flawed, as we now show. 

Let us discuss a definite microscopic model to avoid 
ambiguities. Let H = + JVx^x+i be a Hamilto¬ 

nian of a spin-1/2 chain, where Hx are on-site disorder 
terms with bandwidth VF, Vx^x+i are hopping terms and 


start approaching the value logio[dim(I7) = 6480] ^ 3.5 
of fully delocalized wave functions (cf. the right panel). 
A comparison of histograms at the same values of Ii, but 
with I 2 = 0 (not shown) revealed that the histograms are 
significantly shifted to larger logIPR in the presence of 
the delocalized channel of mobile pairs. 


J is a small coupling satisfying J/W <^1 (large disorder 
condition). 

Let H := [1,..., Li] be a local, finite-size ’’bath”, i.e., 
a part of a chain which is ergodic if isolated from the rest, 
and E := [Li^N] a localized region (the ’’environment”). 
This can arise, e.g., in a disorder realization where all 
Hx ioT 1 < X < Li are equal up to deviations of order 
J, whereas such resonances are absent for x > Li. Now 
split 

H = HbE Hb-e + He (B1) 

where the subscripts indicate the sets of spins on which 
they act (on B only; both on B and its complement E; 
on E only). Following Imbrie [l|, in the localized region 
E we can find a unitary rotation that brings the Hamil¬ 
tonian He into a canonical form of a sum of commuting 
terms, each being a product of ’’logical bit” operators. 
We now apply this rotation to the full Hamiltonian. Only 
the rotation of the terms which involve action on B lead 
to terms that can flip the logical bits. However, their 
amplitude decays as (J/W)^ where R is the distance be¬ 
tween the border of the ergodic region and the site at 
which the logical bit is centered. However, those Hipping 
terms can induce significant hybridization between dif¬ 
ferent values of the logical bits only if their amplitude is 
bigger than the level spacing in the region of diameter 
I/i + i7, which scales as Assuming sufficiently 

strong disorder J/W < 1/2, we define a buffer length Lq 
by 


2-Li-Lo _ 


(B2) 
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The above reasoning implies that logical bits centered 
further than Lq from the edge of the ergodic region are 
not hybridized through their coupling with the ergodic 
region B. To state this in a more formal fashion, we 
can split the chain into B := [1,..., Li + Lq] and E = 
[Li + 1/0 + 1, A/"] and write the (rotated) Hamiltonian as 

H' = + (B3) 

where the terms in acting on both B^E are ei¬ 

ther commuting with (terms diagonal in the logical 
bits), or have norm at most (J/W)^^ (terms originating 
from applying the rotation to Hb-e)- Using (jB2p . this 
allows to conclude that the eigenstates of H' are close to 
products of states in B and eigenstates of the logical bit 
operators in E. In other words, the delocalizing influence 
of the region B is limited to a finite region of length at 
most Lq around H, while many-body localization remains 
robust outside. 

This shows that it is flawed to reason that the ergodic 
region B can be extended to an ever enlarging neighbor¬ 
hood, as the bath would melt its neighboring degrees of 
freedom and incorporate them into the ergodic region. 
In fact, these enlarging regions gradually lose their full 
ergodicity. In particular, random matrix theory cannot 
be applied naively when determining whether further de¬ 
grees of freedom become delocalized as they are coupled 
to the already enlarged region. 

In contrast, our bubble argument in the main text does 
not require such an iterative construction of a growing er¬ 
godic region. The basic difference is that a high-energy 
bubble can exist anywhere in space and is not tied to a 
rare quenched disorder configuration, which remains lo¬ 
calized in space. We thus only need to argue that a hot 
bubble is able to displace itself in a finite series of elemen¬ 
tary steps, thereby remaining ergodic inside. For those 
stages, we do not see reasons to doubt the applicability 
of random matrix-like behavior of matrix elements and 
level spacings. 

It is interesting to note that the above reasoning about 
the robustness of MBL cannot be generalized straight¬ 
forwardly to d > I. The influence of ergodic spots in 
higher dimensional systems thus remains an interesting 
open question. 


Appendix C: Failure of our arguments for single 
particle mobility edges 

It is natural to ask why our argument for delocaliza¬ 
tion by bubbles does not apply to single-particle mod¬ 
els, where we know that mobility edges do exist. In a 
nutshell, the argument does not apply because a single 
particle cannot borrow the energy from other particles 
to become hot and to delocalize at higher energy. Hy¬ 
bridizations between localized and delocalized states (as 
induced by some perturbation to the Hamiltonian) would 
necessarily require to go off-shell, which suppresses the 
coupling between localized states at large distance from 


each other. In contrast, the formation of a mobile bub¬ 
ble in a many-body system does not require the viola¬ 
tion of energy conservation. Let us show that the natu¬ 
ral analog of our argument fails. We have now a single 
particle Hamiltonian H with localized eigenstates \r]) (es¬ 
sentially 177 ) |x)) for some x and extended eigenstates 
|T), separated by a mobility edge at Ec. Upon switching 
on a local perturbation V (say, a change in the random 
potential) we evaluate the possible hybridizations. Of 
course, |? 7 ), \r]') do not hybridize in first order. Also, \r]) 
does typically not hybridize with extended eigenstates 
T. They are separated in energy and the matrix element 
( 77 |U|T) ^ with L the linear size of the system. 

Indeed, this is simply Mott’s argument. (To be precise, 
one should check as well that the sum of all contribu¬ 
tions extended contributions to the change in the state 
1 77 ) does not cause it to delocalize, but we skip this. The 
argument is very similar to what follows.) However, the 
more interesting situation occurs in second order. Let us 
consider a typical term for the mixing of localized states 
77,77' via extended states T: 


E 


I I 

Ejj Erq E^> 


{r]\V\^){'^V\r]') 


(Cl) 


A naive estimate gives that this sum is of order so 
potentially dangerous, but a more careful analysis refutes 
this: let us first change the sum over T by a sum over 77. 
Then, we see that the resulting expression can become 
of order one only if I 77 ) and I 77 ') are close in space and 
their energy difference is comparable to the norm of V. 
This is obviously a case in which we do expect \r]) and 
1 77 ') to hybridize. Hence, we exclude this case and we can 
freely add contributions of localized states to m- Using 
functional calculus, we can recast the resulting expression 
as 


Ef] — Er^> 




Er,-H 


(C2) 


with x{E) a smooth cut-off function (we choose it to be 
a ’Schwarz function’) that vanishes for E < Er^ ^ {Ec — 
Er^)l2 and becomes I for It follows that 


X{H) 

Er,-H 


f{H) 


for a smooth function / of compact support. Now we see 
that f{H) is a local operator in the sense that 

\{x\fi.H)\y)\ < C{k)\x - for any fc > 0 


This is easily proven by Fourier transforming f{H) = 
j f {t), remarking that the Fourier transform of a 

Schwarz function is again Schwarz and using a ’Combes- 
Thomas’ bound to get the locality of on a spatial 

scale proportional to t. For details on the functional anal¬ 
ysis that was used in this argument, see for example jsbj . 
The upshot is that (|C 2 p is given by ^ ^ times an ex¬ 

pression of order ||U|p that decays rapidly in the distance 
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between the localization centers of 17^), 17^'). Hence, in the one-particle theory, the inclusion of extended states in 

perturbation theory does not destroy the localized states. 
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